Delocalization of Wannier-Stark ladders by phonons: tunneling and stretched polarons 
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We study the coherent dynamics of a Holstein polaron in strong electric fields. A detailed 
analytical and numerical analysis shows that even for small hopping constant and weak electron- 
phonon interaction, polaron states can become delocalized if a resonance condition develops between 
the original Wannier-Stark states and the phonon modes, yielding both tunneling and 'stretched' 
polarons. The unusual stretched polarons are characterized by a phonon cloud that trails the 
electron, instead of accompanying it. In general, our novel approach allows us to show that the 
polaron spectrum has a complex nearly-fractal structure, due to the coherent coupling between 
QQ , states in the Cayley tree which describes the relevant Hilbert space. The eigenstates of a finite 

ladder are analyzed in terms of the observable tunneling and optical properties of the system. 
The coupling of an electronic system with a coherent boson field has been a subject of great interest in recent years, 
J> ■ thanks to the availability of high quality materials and intense photon sources. In condensed matter systems, however, 
, a natural intrinsic field of bosonic nature is provided by the lattice vibrations. The coupling to phonons is typically 

■ considered to result in inelastic incoherent scattering for the electron, as the phonons arc difhcult to interrogate 
, separately. This results in an electronic 'open' system which exists immersed in the assumed-incoherent phonon field. 
' Although this description is appropriate for many situations, either because the electron-phonon coupling is weak 

04 , and/or non-resonant, or if the system is at high temperatures, its validity is suspect if the interaction is effectively 

■ strong and the temperature is low. 
Electronic transport in superlattices has yielded a number of interesting phenomena and concepts, including Bloch 

oscillations and Wannier-Stark (WS) ladders 0, as well as negative differential conductance ||] in high DC electric 
^ ', fields, dynamical localization ||], fractional WS ladders under DC and AC electric fields [Q, and resonant magne- 
'"^ topolarons |^] . Similarly, electron-phonon interactions under high electric fields have been investigated in some detail, 
including coupled Bloch-phonon and -plasmon oscillations ^ , as well as phonon-assisted-hopping of an electron on 
O ' a WS ladder A powerful variational treatment of inelastic but coherent quantum transport was presented in 

] while anomalies in transport properties under a resonant condition were studied in and the optical absorption 
K>" ' associated with the resonance of a WS ladder and confined optical phonons was studied in Rcfs. 0,|ll[. 
k>( ] In this paper, we study the effects of coherent coupling of an electron to the phonon system present in its solid at 
; I low temperatures. We are specially interested in the effect of this coupling on the tunneling properties of an electron 
J in an intense electric field region, like the situation achieved in semiconductor superlattices, for example jl^ . Our 
description is non-perturbative for the resonant case with exponentially large Hilbert space, and as such we are able to 
elucidate the role of resonant phonon fields on the otherwise localized electron residing in a WS ladder, for both weak 
and strong coupling. Using a Holstein model for the electron-phonon coupling, and both analytical and numerical 
techniques, we demonstrate that the electron tends to become more extended for increasing coupling with phonons. 
Moreover, we show that resonant coupling results in total delocalization of the polaron for some of the states, and 
the consequent restoring of a 'miniband' structure in the system, even as the applied electric field remains strong. 
Other 'stretched polaron' states, however, are highly degenerate yet strongly localized on one end of a finite structure, 
with a phonon 'cloud' which is nearly detached from the electron (somewhat a precursor of the polaron dissociation 
described in |jl^). We are able to show explicitly that the level spectrum in this regime has a nearly fractal structure 
with rather complex wave functions. On the other hand, when the system is away from the resonant condition, a 
deformed WS ladder results, with electronic wave functions which nevertheless remain localized. We finally analyze 
the tunneling and optical properties one could observe experimentally in such a system. 

We consider a Holstein model fl^ , describing a single electron in a one-dimensional tight-binding lattice which 
interacts locally with dispersionless optical phonons. In addition, we consider that a strong constant electric field is 
applied which generates a WS ladder. The Hamiltonian is then 

Ho = Ej {sjc+cj + t{c+cj+i + c+^iCj)} 
+ ^ Ej a'j'aj + 7 Y.. c+Cj (a+ -f a^- ) , 

where t is the electron hopping constant, oj the phonon frequency, 7 the electron-phonon coupling constant, £j = 
—edEj = —jA the site (or 'well') energy, E the electric field, d the lattice constant, and A the spacing between WS 
'rungs' (we have set 7i = 1). It is known that well- localized WS states will form in the absence of electron-phonon 



interactions, with eigenvalues and eigenfunctions given by Sj, and — Ji-j{2t/edE) \i), where J; is the l-th order 
Bessel function Given the properties of Ji, the WS j-state is clearly localized with a characteristic length 2t/eE 
around site j. It is helpful to introduce the WS creation and annihilation operators, dj — X]i=-oo Ji-j{2t/edE)ci . 
It is easy to show that {df , dj} — Sij, and that the Hamiltonian can be rewritten in terms of dj, dj' . In the case of 
strong electric field (or small hopping constant) {2t/ A ^ 1), the Hamiltonian can be simplified to 



H = Y.. {-Ajd+dj + uia+aj + -/d+dj{a+ + a^) 



A (a. 



)id+d,+,+d+d,)} , 



(2) 



where the effective hopping becomes A = '^t/A. Notice that phonon-assisted hopping develops between the WS states 
via the last term. Eq. is similar to that in but here, as the electron jumps it can not only emit (or absorb) 
a phonon on its current location, but also on the neighboring site. This more natural description results in a much 
larger relevant Hilbert space for this problem. 

Consider the process of an electron jumping between WS states and creating or annihilating a phonon with 
amplitude A, as per the last term in (^. The relevant Hilbert space can be described for n = 3 via the basis: 
|1;000),|2;010),|2;100),|3;011),|3;101),|3;020), and |3;110). Each state |j; mi, ma, ...) is indexed by the electron 
position J, and ruk refers to the number of phonons on site k. This construction yields a Hilbert space of dimension 
iV — 1 = 2" — 1, where n is the total number of sites in the chain. It is interesting to note that the structure of this 
vector space is that of a Cay ley tree where each j-state links to two at + 1, in general with non-symmetric weights 
in all branches (isj . A similar structure of the vector space was found in |^ . 

Under the resonance condition, i.e., A = u, all the basis-states in the Hilbert space described above have the same 
energy, and one would expect that the off-diagonal matrix elements in (^) would break the degeneracy. One would 
expect that a 'miniband' of extended states might form, whenever the resonance condition is reached. For this case, 
Hamiltonian (||) in the above basis of n = 3 takes the form 
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with ever larger blocks for larger n values. Notice that there are elements with somewhat unusual values, — v^A, 
associated with the higher number of phonons in state |3; 020). A similar structure is seen for all values of n, and as 
these terms appear only sporadically, they have only a small influence on the overall level structure. In particular, one 
expects that there would be little change in the physics in the limit of large N, if we replace — \/2A terms with —A. 
This assumption is in fact confirmed explicitly by numerical calculations. Moreover, this 'symmetrized' Hamiltonian 
Hsym allows one to better analyze the structure of the solutions oi H, as we show below, while exhibiting some of the 
same generic features and fully self-similar properties. 

Symmetrized Hamiltonian. By using the block decomposition of the determinant, det 



A B 
C D 



det{D) det{A- 



•4-1 = 0, 



BD^^C) , one can find that the eigenvalues of Hsym are determined by the equation Sg^'^e^^'^e^^^ ■ 
where Eq = e, e is the energy eigenvalue, and e^+i = e — 2}? jsh ■ Here we have made a constant energy shift of A. 
One can then write ei^ as a continuous fraction in fc steps 



2A^ 



2A^ 



(4) 



and the eigenvalues of H^ym are given by [|TT[ 

£fc,, =2V2Acos(^^) , 



(5) 



where = 1, 2, n\l = \, k. The limiting bandwidth of the spectrum is 4-\/2A. 

Using Hsym, one can arrive at an estimate for the degeneracies for the first few eigenvalues in the limit of large N , 
given by the equation e'^^/^(£^ — 2A^)^/*(£^ — 4A^)^/^ = 0. One then has eigenvalues e = with degeneracy 3iV/8, 



£ = ±a/2A with degeneracy and e — ±2A with degeneracy -/V/8, in agreement with our expUcit calculations 

for finite N . Notice that increasing lattice or chain size produces both new energy eigenvalues within the limiting 
band and additional states which augment the degeneracy of some of the old eigenvalues. The overall structure of the 
spectrum remains the same and all the eigenvalues of the shorter chain are also in the spectrum of the longer one [0 . 

Figure 1 shows typical numerical results for the eigenvalues of the physical Hamiltonian H in Eq. while the 
inset shows the results for the symmetrized version Hsym (both for n = 8) . Our analytical results match exactly those 
shown in the inset, which in turn mimic well the results for H in the main panel, although small differences appear 
due to the sporadic — v^A terms. The self-similar structure of the spectrum in the inset is the manifestation of the 
full self-similarity of the Hamiltonian Hsym (while it is only approximate for H). One can see that the original highly 
degenerate Hilbert basis is split into a semi-continuous 'miniband' by the off-diagonal elements in H (and which for 
Hsym yields Eq. ^ for large N). Notice, however, that there are large residual degeneracies for the eigenvalues of H 
in the center of the miniband. The right panel in Fig. 1 shows the density of states of the H system, exhibiting also 
an approximate self-similarity. 

The 'electronic wave functions' for various sample eigenstates are shown in Fig. 2, P^{j) — X]{m} l^j {m}!^' "^here 
the e-eigenstate is given in terms of the basis above by = ^^^j {r?x}b ' {™})- C'ne can observe that the states 
in regions of high degeneracy, e = and e/A w —1.6 (where states (A) and (C) are indicated), are localized around 
the right end of the structure, j w 8. However, non-degenerate states in the band, such as the one labeled (B), result 
to be fully extended over the structure. 

Notice that in this inherently coupled electron-phonon system, the calculated eigenstates are a coherent mixture of 
electron and phonons on different sites. The amplitudes in Fig. 2 have been projected over all phonon components. 
Similarly, in order to illustrate the physics of these states better, one can exhibit the average phonon number for 
various states in Fig. 3. One can see clearly that all states at the degenerate center of the band have many more 
phonons (~ 7) than those non-degenerate states residing higher in the band. It appears, in fact, that when the 
electron interacts with many phonons, it localizes, as one would perhaps intuitively expect for polaron self-trapping 
states, while states with fewer phonons are less localized and 'propagate' better. One can provide an estimate for the 
average number of phonons for extended states by the following argument: The average phonon number is given by 
{N) — PiNi, where Pi is the probability for an electron to be at site i, and Ni is the number of phonons in this 
state. For an extended state. Pi ~ where n is the number of sites. A state with an electron at site i can be 
obtained by the electron hopping i — 1 steps from the first site, while emitting Ni = i — \ phonons in the process. 
We then obtain that the average number of phonons for an extended states is {N) ~ X]i(* — l)/?^=('^■— l)/2. This 
result matches well the numerical results, where the most extended states have {N) ~ 3.5 in Fig. 3, such as the state 
labeled (B). 

One can also inspect the spatial distribution of phonons for each state, as in the inset of Fig. 2, which shows the 
corresponding phonon distribution for several eigenstates, N^{j) = j^^,^ {m-}\'^''^h indicating how the phonons 
distribute over the chain. The degenerate states, such as (A) and (C), with the electron localized to the right of the 
structure, have nearly one phonon per site throughout the chain, indicating that these polarons have a 'stretched' 
phonon cloud away from the charge, seemingly in anticipation of the total polaron dissociation that other authors 
have described for extended chains under high electric fields |Q. In contrast, the most delocalized states (such as 
(B) in Fig. 2) show a phonon amplitude much less than unity throughout, yielding an overall low phonon count, and 
a cloud that in general 'accompanies' the charge. 

Interband optical absorption. The spectrum of optical absorption associated with transitions between valence and 
conduction bands provides information about the density of states and the states' charge distribution. The absorption 
intensity K is calculated from K = '^\Wv\'^5{£^ + Eg — tu), where Eg is the valence and conduction band gap, i/ 
refers to the various polaron states, and Wi, is the transition matrix between valence band and conduction band 
states. Under the strong field present in the system, the hole is localized the most because of its large effective mass. 
Thus we ignore its negligible tunneling and coupling to phonons and assume it localized at site 1. Correspondingly, 
there is conservation of phonon number and jW^p cx Pu{^) [@- In Fig. 4(a) we show the absorption intensity vs. 
phonon frequency. Once again, we see that the non-degenerate states, e/A ~ ±0.6, have the largest amplitude at site 
1 and will then contribute the most to the absorption spectrum. In contrast, all highly degenerate states at e = are 
strongly localized at the right end of the chain and make a vanishing contribution to K. Such strong modulation of 
K would be easily accessible to experiments. 

Transport properties. As described above, tunneling experiments under strong electric fields are possible in semicon- 
ductor superlattices and other systems, and it is of interest to investigate the role of the resonance condition. The tran- 
sition probability P can be calculated in a S'-matrix formalism. P= |rpandr= {e f, R\S\£i, L) . Here i? and i refer to 
the right and left leads. In the wide band limit, T cx F / dtidt2e''^^f*^~^^*'^n\G"{t2-ti)\l) = E,. r(ri|i/)(i^|l)(5(£/-ej). 



where F describes the electron interaction with contacts, and is the retarded Green function connecting both ends 
of the structure (the site j = I and j ~ n) We use D{i') = |(n|z/)(z/|l)p = Py(l)Py(n) to describe the relative 
contribution of various eigenstates v to the electronic transport, as the tunneling probability is proportional to the 
density of states in the leads and D{i^) |19|. Figure 4(b) shows the quantity D{e). We can see that all the states 
at the center of the band contribute zero to the transport amplitude through the chain, and that this behavior is 
exhibited by all the high peaks in the DOS, clearly consistent with their spatially localized-charge nature. On the 
other hand, D shows large values for extended states, confirming in fact that the non-degenerate states in the mini- 
band are extended. The variations shown in D would then be reflected in strong amplitude modulations within each 
of the phonon replicas in tunneling experiments ]20|] , whenever the resonance condition is reached. It is also clear 
that despite the fact that K and D correspond to completely different physical quantities, their energy behavior is 
surprisingly similar. This arises because of the large electric field asymmetry, resulting in no states with mostly/only 
large electronic amplitude in the first site (i.e., localized towards the left of the chain). Once a state is extended and 
has non-vanishing amplitudes at site 1, it also does at site n. 

So far we have mainly described the resonant case. When the system is away from the resonant condition, a deformed 
lattice with substructure is formed. We are also able to study the system beyond small hopping and weak coupling 
regime by mapping it via the Lang-Firsov canonical transformation to the case we have studied, with renormalized 
hopping and coupling constants. In this approach, we can provide additional physical insights into some interesting 
results found in We will report the details elsewhere ]l7t . 

Although our model is ID | p2[ , we expect the qualitative features and main conclusions to be valid for semiconductor 
superlattices (SLs) where electron-optical-phonon coupling is highly anisotropic. In that case, it would be possible to 
neglect the in-plane scattering of electrons by phonons, while the strong interaction assists in the electron hopping. 
In fact, for SLs with low in-plane disorder, the transitions between different layers involve states near the bottom 
of the 2D electron subbands and phonons with small in-plane momenta. Hence, the system would exhibit quasi- ID 
behavior, and in-plane scatterings are a small correction. Similarly, we could adapt our model to the description of 
polarons in polymer chains, such as those described in Ref. [|l3|, although there the electron- phonon interaction is of 
different strength than in solid SLs, and tunneling experiments are rather more challenging. 

We have studied the coherent dynamics of Holstein polarons in a strong electric field. It is found that with the help 
of phonons, each WS rung develops into a 'miniband' under the resonance condition. This miniband shows a nearly 
self-similar structure, which is inherited from the full self-similarity of the Hamiltonian Hsym for a symmetric Cayley 
tree. Although the phonons can help the electron jump from one WS state to another, the phonons can also prevent 
the electron from propagating. If too many phonons are involved, this results in highly degenerate states with the 
electron localized at one end of the structure while the phonons are pulled away from it (the 'stretched' polaron). The 
miniband structure and its modulations in density of states and other characteristics are manifested in both transport 
and optical properties of the system, which we anticipate can be observed in experiments at low temperatures |12|. 
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FIG. 1. Energy spectrum for H 

structure of spectrum for Hsym, a 
A to C relate to Fig. 2. Here, A = 



n = 8. Notice clear fractal 
metric about e = 0. Labels 
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FIG. 2. Electronic amplitudes P{j) at each site j for several states, A to C, as indicated in Fig. 1. Notice A and C are 
localized near j ~ 8. Panel shows phonon content N{j) for same states. A and C have large phonon counts throughout and 
away from the electron, resulting in 'stretched' polarons. 









I 










r 




(A) 




n 


) 




(C) 








\ 




I 










(B) 







■i H 1 1 1 1 1 1 1 1 1 1 1 1 

50 100 150 200 250 300 

State index 

FIG. 3. Total phonon number for each eigenstate. Highest values plateau for the most localized electronic amplitudes, A 
and C, while lowest counts correspond to extended electronic states, such as B. 
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FIG. 4. (a) Optical absorption spectrum K within miniband with level broadening added, (b) Tunneling amplitude D. Both 
cases for spectrum in Fig. 1. 



